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1 J  INTRODUCTION 

v) ihe  Air  Force  and  other  Department  of  Defense  agencies  have  for  many  years 
supported  the  development  and  improvement  of  infrared  radiance  and  transmis¬ 
sion  models  of  the  atmosphere.  These  models  are  designed  to  simulate  the  ra¬ 
diative  properties  of  the  atmosphere  to  provide  predictions  for  Air  Force  and 
other  Department  of  Defense  system  design  and  surveillance  programs.  Compari¬ 
son  with  available  experimental  measurements  provides  greater  understanding  of 
the  atmosphere  and  serves  to  check  the  adequacy  of  the  models. 

This  report  is  intended  to  provide  a  User's  Guide  to  such  a  radiance 
model,  the  AFGL/Visidyne  High  Altitude  Infrared  Radiance  Model.  This  model  is 
designed  to  partially  simulate  the  earth's  atmospheric  radiance  in  an  earth's 
limb  viewing  geometry.  The  nominal  spectral  interval  included  at  present  lies 
between  2.0  and  40  micrometers.  Emphasis  is  placed  on  radiation  originating 
between  60  and  250  km,  although  the  model  has  been  used  to  provide  approximate 
radiance  values  for  tangent  heights  as  low  as  50  km  and  as  high  as  500  km. 

Various  versions  of  this  model  have  been  described  by  Degges^* and  by 
Degges  and  Smith  Degges^  gave  a  listing  of  the  program  with  a  very 
brief  description  of  program  logical  flow.  This  model  was  restricted  to  alti¬ 
tudes  above  60  km  altitude,  and  included  the  infrared  active  gasses,  methane, 
carbon  dioxide,  water  vapor,  nitrous  oxide,  nitric  oxide,  and  ozone.  The  spec¬ 
tral  region  covered  was  between  five  and  twenty-five  micrometers,  and  only  a 
few  vibrational  levels  were  included  for  each  species.  At  the  time  of  this 
report,  there  were  very  few  experimental  observations  of  the  abundance  of  in¬ 
frared  emmitting  molecules  at  the  altitudes  required  by  the  program.  This  work 
included  an  additional  computer  program  to  estimate  diurnal  variation  in  the 
chemical  species  of  interest,  to  investigate  the  effects  of  changes  in  assumed 
eddy  mixing  coefficients,  and  to  estimate  the  effects  of  new  determinations  of 
chemical  rate  coefficients.  This  work  on  the  chemistry  has  not  been  updated 
since,  and  we  currently  prefer  to  rely  on  more  up-to-date  published  calcula¬ 
tions  and  observational  determinations  of  chemical  species  abundances. 

The  principal  areas  of  research  and  program  revision  reported  by  Degges^ 
and  Degges  and  Smith  ^  included  a  further  study  of  band  radiance  modeling  and 
the  effects  of  the  line  shape  on  radiance  computations.  In  addition  to  some 
exploratory  calculations  using  a  Voigt  line  profile  instead  of  a  Doppler  line 
shape,  some  asymptotic  analytical  results  were  obtained  which  provided  an  es¬ 
timate  of  the  accuracy  of  the  numerically  obtained  band  functions  at  large  op- 
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tical  thicknesses.  It  was  found  that  when  a  Voigt  profile  is  used,  the  connec¬ 
tion  between  the  band  models  and  functions  for  a  single  line  is  more  readily 
seen  than  when  a  Doppler  line  shape  is  used.  The  program  described  here  con¬ 
tinues  to  be  based  on  a  model  with  Doppler  line  shape.  This  restricts  the 
lower  altitude  for  accurate  quantitative  results  to  about  60  km. 

Since  that  work,  there  have  been  several  modifications  to  the  computer  pro¬ 
gram.  Additional  carbon  dioxide  and  ozone  bands  have  been  added,  the  option 
of  computing  radiative  cooling  (or  heating)  rates  has  been  added,  and  some  as¬ 
pects  of  the  computational  method  used  in  computing  carbon  dioxide  and  ozone 
vibrational  populations  have  been  revised.  The  addition  of  the  carbon  dioxide 
bands  was  straight-forward,  but  more  effort  was  required  in  adding  ozone  bands 
to  the  ten  to  twelve  micrometer  "window"  region. 

In  addition  to  model  changes,  we  have  used  the  U.S.  Standard  Atmosphere 
f  41 

Supplements,  1966 1  J  to  provide  a  set  of  atmospheric  input  data  which  can  be 
used  to  demonstrate  some  of  the  large  scale  effects  of  temporal  2nd  latitudi¬ 
nal  variations  of  the  atmopshere.  Section  2  presents  this  atmospheric  data 
set.  Section  3  provides  a  description  of  use  of  the  computer  implementation 
of  the  model  on  the  AFGL  computer  system.  Section  4  describes  the  effort  to 
expand  the  ozone  radiance  part  of  the  model. 
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2.  DATA  FOR  MOLECULAR  CONCENTRATIONS  AND  INFRARED  RADIANCES 

The  principal  data  required  for  calculation  of  non-equilibrium  infrared 
radiances  include  the  concentrations  of  atmospheric  species ,  the  rate  coeffi¬ 
cients  for  collisional  excitation  of  infrared  emitting  vibrational  levels  and 
infrared  band  strengths  and  solar  flux  values.  In  previous  work  we  have  usually 
used  one  of  the  U.S.  Standard  Atmospheres  (1962^  or  1976^)  for  information 
on  molecular  nitrogen  and  oxygen  and  atomic  oxygen,  and  other  sources  for  the 
concentration  of  the  infrared  emitting  species,  we  follow  this  practice  here, 
in  providing  a  number  of  data  sets  for  the  principal  atmospheric  species  and 
what  is  essentially  a  single  set  of  concentrations  for  002,  H20,  NO,  0^,  and  CO. 
We  retain  the  use  of  concentrations  based  on  the  U.S.  Standard  Atmosphere, 
1976*^,  and  supplement  it  with  24  other  data  sets  based  on  the  U.S.  Standard 

r  41 

Atmosphere  Supplements,  1966 1  . 

f  41 

The  U.S.  Standard  Atmosphere  Supplements,  1966 1  J  provides  model  atmospheres 
for  altitudes  up  to  120  km  at  15  degree  latitude  intervals  between  15  and  60 
degrees  north  latitude.  The  tropic  model  at  15  degrees  latitude  is  assumed 
constant  year  round.  Winter  and  summer  models  are  provided  for  the  other 
latitudes,  and  a  spring  or  fall  model  is  added  at  45  degrees.  High  altitude 
atmospheres  are  given  for  combinations  of  winter,  summer,  or  spring/fall 
conditions  and  exospheric  temperatures  between  500  and  2500  K  to  simulate 
variations  in  the  solar  cycle.  We  have  selected  a  set  of  these  upper  atmosphere 
models  (600,  1000,  1500  K  exospheric  temperatures)  to  make  possible  calculations 
that  can  illustrate  aspects  of  atmospheric  variability.  We  use  temperature 
profiles  and  sea  level  pressure  to  compute  concentrations  of  N2  and  02  up  to  90 
or  110  km,  assuming  constant  mixing  up  to  110  km  and  match  the  high  altitude 
profiles  with  the  computed  lower  atmosphere  concentrations. 

Table  1  lists  temperature  profiles  below  90  km  for  eight  atmospheres  from 
the  1966  Supplements.  Table  2  lists  nine  high  altitude  temperature  and 
concentration  profiles,  again  from  the  1966  Supplements.  Table  3  gives  the 
temperature  profiles  for  the  24  atmospheres  produced  by  combining  the 
information  in  Tables  1  and  2,  as  well  as  the  1976  Standard  Temperature 
Profile^.  Table  4  gives  day  and  night  atomic  oxygen  profiles  for  the  nine 
atmospheres  of  Table  2.  These  are  the  same  as  those  used  with  the  1976  Standard 
Atmosphere  up  to  120  km.  Table  5  lists  high  altitude  concentrations  for  the 
1976  Standard  Atmosphere. 
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For  the  other  species,  we  use  the  model  of  Degges  and  Smith ^  ^ ^  except  for 
water  vapor  and  carbon  monoxide  mixing  ratios  from  the  calculations  by  Allen 
et  al.^.  This  model  assumes  a  constant  mixing  rate  for  COj  up  to  110  km, 
with  diffusive  equilibrium  above  that  altitude.  This  results  in  daytime 

concentrations  that  are  too  large,  probably  by  a  factor  of  two  or  more  above 

150  km.  Our  nitric  oxide  concentrations  are  based  on  the  twilight 
measurements  of  Baker  et  al.^  and  on  the  calculations  of  Oran  et  al. 

The  "Day  Time"  nitric  oxide  concentrations  are  better  described  as  referring 
to  high  latitude  conditions.  The  "Night  Time"  ones  are  closer  to  mid-latitude 
conditions. 

Below  60  km,  we  use  the  ozone  model  in  the  U.  S.  Standard  Atmosphere, 

1976.  Above  that  altitude,  our  models  rely  heavily  on  measurements  by 

Llewellyn  and  Witt^®^  and  Hays  and  Roble^*^.  The  calculations  of  Thomas  and 
f  121 

Bowman1  have  been  used  as  a  guide  to  predict  diurnal  variations  in  abun¬ 

dances.  Diffusive  equilibrium  is  assumed  above  100  km. 

These  concentration  profiles  should  be  viewed  as  models,  and  not 
considered  to  represent  current  knowledge  of  the  atmosphere.  Work  is  in 
progress  at  AFGL  to  provide  better  concentration  profiles  based  on  the 
Mesospheric  Explorer  Satellite  data.  The  current  models  will  then  be  revised 
to  reflect  this  new  experimental  information. 

Some  rates  for  collisional  de-excitation  of  vibrational  levels  have  been 
taken  from  Taylor We  continue  to  use  the  nitric  oxide  quenching  by 
atomic  oxygen  rates  given  by  Glanzer  and  Troe^^.  The  water  vapor  calcula¬ 
tions  have  been  revised  to  use  the  quenching  rates  of  Hovis  and  Moore ^  ^  and 
Finzi  et  alJ17^.  We  use  new  rates  suggested  by  R.  SharnaJ1®^  for  OC>2  15  /m. 
bands. 

Tables  6  through  10  list  the  molecular  bands  and  related  physical  data  now 

included  in  the  programs.  Band  strengths  for  water  vapor,  carbon  dioxide, 

ozone,  and  carbon  monoxide  are  those  used  in  the  AFGL  absorption  line  atlas 

r  19] 

first  described  by  McClatchey  et  al.1  .  Band  strengths  for  nitric  oxide  are 
approximately  those  of  Abels  and  Shaw^2^  and  Schurin  and  Ellis^2*^.  Infrared 
solar  fluxes  at  band  centers  are  interpolated  from  the  tabulation  of 
Johnson^ 22 ^ . 


TABLE  1.  DEFININ';  ALTITUDES,  TTMPERATOKES , AND  SEA  LE\/EL  PRESSURE  FOR  ITWER 
ALTITUDE  MODEL  ATMOSPHERES. 


TEMP  45  JANURY  004  U.S.  STND.  ATM.  1966 

0.00  272.15  3.01  261.65  10.02  219.65  19.06  215.15 
27.11  215.15  32.16  219.15  47.35  265.65  52.43  265.65 
64.65  241.65  85.12  199.65  91.29  199.65  101.60  227.65 


TABLE  1  DEFINING  ALTITUDES,  TEMPERATURES, AND  SEA  LEVEL  PRESSURE  FOR  LOWER  ALTITUDE 
MODEL  ATMOSPHERES. 


PRES  60  001Y  008 

10l0'°0 


TABLE  2.  TEMPERATURE  AND  MOLECULAR  NUMBER  DENSITIES  OF  HIGH 
ALTITUDE  MODEL  ATMOSPHERES. 
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TABLE  5.  HIGH  ALTITUDE  MOLECULAR  NUMBER  DENSITIES  FOR  U.S.  STANDARD 


ATMOSPHERE,  1976 


U.  S.  STANDARD  ATMOSPHERE.  1976 


86.00  1 . 1 300E+ 1 4 

88.00  7 . 91 50E+1 3 

9o.oo  5.5470E+13 

93.00  3 . 806OE+1 3 

94.00  2 ■ 71 50E+1 3 

96.00  1.8940E+13 

98.00  1.3200E+13 

100.00  9 . 2 1 OOE+1 2 

102.00  6.50BOE+12 

104.00  4 . 6090E+1 2 

106.00  3 . 27  30E  +  1 2 

108.00  2 . 32706+12 

110.00  1 . 64 1 0E+ 1 2 

112.00  1 . 1 580E+1 2 

114.00  8.4220E+11 

116.00  6 . 2850E+ 1 1 

118. 00  4 . 7940E+1 1 

120.00  3 . 7260E+ 1 1 

122.00  2.9470E+11 

124.00  2 . 3680E+1 1 

126.00  1 . 9300E+1 1 

128.00  1.59206+11 

130.00  1 . 32606+1 1 

132.00  1 .11606+11 

134.00  9.46006+10 

136.00  8 . 0800E+ 1 0 

138.00  6 . 94706+10 

140.00  6.00906+10 

142.00  5.22506+10 

144.00  4.5CS0E+10 

146.00  4 . 00706+1 0 

148.00  3.53106+10 

150.00  3 . 1 2406+1 0 

160.00  1 . 77406+ 1 0 

170.00  1 . 0700E+ 1 0 

180.00  6 . 74006+09 

130.00  4.3850E+09 

200.00  2.92501+0') 

220.00  1.37306+09 

240.00  6.7780E+08 

260.00  3 . 4590E+08 

280.00  1 . 806QE+0Q 

300.00  9 . 59306+0? 

350.00  2.06906+07 

400,00  4 . 6690E+0G 

450.00  1 . 0860E+06 

500.00  2 . 59206+05 

550.00  6.32306+04 

600.00  1 . 57506+04 

650.00  4 . 00306+03 

700.00  1.0380E+03 

800.00  7.37706+01 

900.00  5.641 06  +  00 

1000.00  4.6260E-01 

# 


3.03106+13  8. 60006*1 0 

2.  1  2006+13  1.5U0E+11 

1 . 4?90E+ 1 3  2.4430E+11 

1 . 0270E+1 3  3.4340E+11 

7 . 06006+1 2  4.1590E+11 

4.8010E+12  4.4710E+11 

3 . 2300E  +  1 2  4.4760F.+  U 

2.151 OE  +  1 2  4 . 2980E+1  1 

1.4300E+12  4 . 0070E+1 1 

9.4340E+11  3.6190E-M1 

6. 1 890E+1 t  3. 1880E+1 1 
4 . 0450E+1 1  2.74806+11 

2.62106+11  2.3030E+11 

1.7060E+11  1 . 88906+1 1 

1 . 1 560E+1 1  1 . 5650E+1 1 

8 . 1 200E+1 0  1 . 3050E+1 1 

5 . 8920E  +  1 0  1 . 0960E* 1 1 

4 . 39506+ 1 0  9 . 27506+1 0 

3 . 3600E+1 0  7 . 9250E+ 1 0 

2 . 6250E  +  1 0  6.8400E+10 

2.0870E+10  5.956OE+10 

1 . 6830E+ 1 0  5.2290E+10 

1 . 3750E+ 1 0  4 . 6250E+ 1 0 

1 . 1 340E  + 1 0  4.11 0OE  +  1 0 

9.4440E+09  3 . 6980E+1 0 

7 . 9270E+09  3.32006+10 

6.70206+09  3.00406+10 

5 . 70206+09  2 . 72906+ 1 0 

4.881 0E+09  2 . 48906+1 0 

4 . 1 990E+09  2.2V80E+10 

3.6310E+09  2.0920E+10 

3 . 1 530E+09  1 . 9270E+ 1 0 

2 . 7500E+09  1  .  7800E  +  1 0 

1.46006+09  1  .  23806+1  0 

8.2770E+08  8.9960E+09 

4.921  OE+OB  6.7470E+09 
3,031 OE+OB  5 . 1 3 1 0E+09 
1.91 806+08  4.05006+09 

8 . 1 450E+07  2.5730E+09 

3.6530E+07  1 .69506+09 

1.7000E+07  1  . 1 4  30E  +  09 

8.11 00E+06  7.03406+08 

3.94206+06  5.4330E+08 

8.85006+05  2.24906+08 

1.25206+05  9.56406+07 

2.36606+04  4.1640E+07 

4.60706+03  1.63606+07 

9.19606+02  8.2000E+06 

1.8600E+02  3 . 7070E+06 

3.93206+01  1.6950E+06 

8.41006+00  7.8400E+05 

4. 10506-01  1.73206+05 

2.17706-02  3.98906+04 

1.25106-03  9.5620E+03 


TABLE  6:  WATER  VAPOR  BAND  DATA 


VIBRATIONAL 

TRANSITION 

BAND  CENTER 

Upper 

STATE 

LOWER  ' 
STATE 

WAVELENGTH 

(um) 

010 

000 

6.27 

020 

000 

3.17 

100 

000 

2.73 

001 

000 

2.66 

020 

010 

6.42 

100 

010 

4.85 

001 

010 

4.63 

030 

010 

3.26 

no 

010 

2.75 

Oil 

010 

2.68 

030 

020 

6.60 

030 

000 

2.14 

no 

000 

1.91 

Oil 

000 

1.88 

BAND  STRENGTH 
AT  296 °K 

ip  i 

(mol  cm  cm  ) 

SOLAR  FLUX  AT 
BAND  CENTER 

-2  1 

(photons  cm  um 

1 . 06 (-17)* 

5.60(15) 

6 . 58 ( -20 ) 

3.50(16) 

3.62( -19) 

5.10(16) 

7 .99( -18) 

5.40(16) 

8.6K-21) 

5.70(15) 

8.90(-23) 

1.20(16) 

4 . 90 ( -22) 

1.50(16) 

7 . 99 ( -23 ) 

3.30(16) 

1  -  50 ( -22) 

5.00(16) 

3 . 30( -21 ) 

5.40(16) 

7 .00( -24) 

4.90(15) 

2.00(-22) 

9.30(16) 

1 . 83 ( -20) 

1.20(17) 

9.16M9) 

1.25(17) 

*Number  in  parentheses  is  power  of  ten. 
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TABLE  7:  CARBON  DIOXIDE  BAND  DATA 


VIBRATIONAL  TRANSITION 

BAND  CENTER 
WAVELENGTH 
(Mm) 

- UPPER 

STATE 

LOWER 

STATE 

01101 

00001 

14.98 

10002 

01101 

16.18 

02201 

01101 

14.98 

10001 

01101 

13.87 

11102 

00001 

5.17 

11102 

10002 

15.45 

11102 

02201 

16.74 

11102 

10001 

18.37 

03301 

02201 

14.97 

11101 

00001 

4.81 

11101 

10002 

12.64 

11101 

02201 

13.48 

11101 

10001 

14.52 

00011 

0000 1 

4.26 

00011 

10002 

9.40 

00011 

10001 

10.41 

01111 

01101 

4.28 

10012 

00001 

2.77 

10012 

10002 

4.30 

02211 

02201 

4.30 

10011 

00001 

2.69 

10011 

10001 

4.30 

20013 

00001 

2.06 

20013 

10002 

2.80 

20012 

00001 

2.01 

20012 

10002 

2.71 

20012 

10001 

2.79 

20011 

00001 

1.96 

20011 

10001 

2.69 

BAND  STRENGTH  SOLAR  FLUX  AT 


AT  296*K 

1  2  1 
(mol'  cnrcnf  ) 

BAND  CENTER 

-2  -1 

(photons  cm'  pm"  ) 

8.26( -18)* 

4.60(14) 

1 . 44 ( - 19 ) 

3.70(14) 

6.49(-19) 

4.60(14) 

1 .85 ( -19) 

5.80(14) 

1 . 06 ( -21 ) 

9.45(15) 

2 .22 ( -20) 

4.20(14) 

5.2K-21) 

3.35(14) 

? . 72 ( -22 ) 

2.55(14) 

3.82( -20) 

4.60(14) 

4 . 87 ( -21 ) 

1.15(15) 

1 . 12 ( -21 ) 

7.00(14) 

7 . 90 ( -21 ) 

6.30(14) 

1 . 49 ( -20) 

5.00(14) 

9 . 60 (-17) 

1.60(16) 

9 . 03 ( -22 ) 

1.78(15) 

6 . 87 ( -22 ) 

1.31(15) 

7 . 66 ( -18) 

1.60(16) 

1 .04 ( -18 ) 

4.95(16) 

1 . 93 ( -19) 

1.60(16) 

3 . 08 ( - 19 ) 

1.60(16) 

1 . 50( -18) 

5.35(16) 

1 . 18 ( -19 ) 

1.60(16) 

8 . 07 ( -21 ) 

1.02(17) 

3 . 38 (-21 ) 

4.80(16) 

3 . 50 ( -20) 

1.07(17) 

4 . 24 ( -21 ) 

5.20(16) 

1.78(-21) 

4.90(16) 

1 . 12 ( -20) 

1.14(17) 

3 . 50 (-21) 

5.30(16) 

*Number  in  parentheses  is  power  of  ten 
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TABLE  8:  NITRIC  OXIDE  BAND  DATA 


BAND  STRENGTH  SOLAR  FLUX  AT 
VIBRATIONAL  TRANSITION  BAND  CENTER  AT  296#K  BAND  CENTER 

- UPPER  LOWER”  WAVELENGTH  \  \  21 

STATE  STATE  (wm)  (molcnrcm" i)  (photons  cm"  urn-1) 


TABLE  9:  OZONE  BAND  DATA 


VIBRATIONAL  TRANSITION 


STATE 


STATE 


BAND  CENTER 
WAVELENGTH 
(pm) 


BAND  STRENGTH 
AT  296#K 

-12  1 
no!  cm  cm  ) 


SOLAR  FLUX  AT 
BAND  CENTER 

(photons  cm  pm  ) 


010 

000 

14.27 

6.28 

X 

10  17 

10  9 

10_20 

in 

-20 

10-19 

10-20 

10-20 

10-19 

10-2l 

10  £L 

5.30  x 

001 

000 

9.60 

1.29 

X 

1.68  x 

100 

000 

9.07 

6.71 

X 

1.97  x 

020 

000 

14.32 

4.16 

X 

5.20  x 

Oil 

000 

5.79 

5.37 

X 

7.50  x 

Oil 

010 

9.75 

4.50 

X 

1.60  x 

110 

000 

5.57 

2.27 

X 

7.70  x 

no 

010 

9.13 

1.10 

X 

2.25  x 

002 

000 

4.86 

1.11 

X 

1.12  x 

002 

001 

9.84 

1.74 

X 

1.55  x 

101 

000 

4.74 

1.13 

X 

1.20  x 

101 

001 

9.36 

3.22 

X 

1.80  x 

101 

100 

9.92 

6.15 

X 

10-20 

O'20 

8-2! 

-20 

10-20 

10-22 

10-22 

10 

-22 

10  cc 

O"21 

o-19 

-21 

10-21 

10-23 

10 

-25 

10  Yi 

1.50  x 

200 

000 

4.54 

3.00 

X 

1.35  x 

200 

100 

9.11 

6.30 

X 

1.92  x 

111 

000 

3.59 

2.32 

X 

2.55  x 

111 

010 

4.80 

3.85 

X 

1.16  x 

111 

001 

5.74 

1.52 

X 

7.10  x 

111 

100 

5.94 

2.65 

X 

6.50  x 

111 

Oil 

9.45 

1.65 

X 

1.75  x 

111 

no 

10.11 

2.50 

X 

1.43  x 

003 

000 

3.29 

1.11 

X 

3.20  x 

003 

001 

5.00 

2.23 

X 

1.04  x 

003 

002 

10.17 

1.96 

X 

1.41  x 

004 

003 

10.50 

2.12 

X 

1.38  x 

005 

004 

10.85 

2.66 

X 

1.38  x 

006 

005 

11.20 

3.61 

X 

10-27 

"1.-29 

10  ii 

1.38  x 

007 

006 

11.60 

5.48 

X 

1.38  x 

008 

007 

12.02 

9.50 

X 

10-32 

10  i£ 

1.38  x 

009 

008 

12.47 

1.87 

X 

1.38  x 
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TABLE  10:  CARBON  MONOXIDE  BAND  DATA 

BAND  STRENGTH  SOLAR  FLUX  AT 


VIBRATIONAL 

TRANSITION 

BAND  CENTER 

UPPER- 

LOWER 

WAVELENGTH 

STATE 

STATE 

(Mm) 

1 

0 

4.67 

2 

0 

2.35 

2 

1 

4.72 

AT  296  K  BAND  CENTER 

12-1  -2  -1 
(mol"  cm  cm"  )  (photons  cm  um~  ) 

9.70  x  10'i!  1.25  x  lO.1® 

6.99  x  10",,  7.40  x  10}J 

5.80  x  10  1.20  x  10lb 
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3.  USING  THE  COMPUTER  MODEL 

Ihe  nigh  altitude  radiance  program  oonsists  of  four  program  sets.  The  first 
includes  five  programs  that  compute  vibrational  populations,  BGND2  for  CO 2, 
bGND3  for  H^O,  BNGD4  for  NO,  ENGD6  for  0^,  and  BNGD7  for  CO.  The  other  three 
are  BGND9,  SPCTRA  and  TRYPEN.  Each  program  may  be  run  separately  or  in 
appropriate  combination  with  one  or  more  of  the  others.  A  general  flow  chart 
for  tne  high  altitude  radiance  program  is  given  in  Figure  1,  which  shows  the 
basic  routines  of  the  oomputer  program.  We  use  the  water  vapor  program  BGND3  to 
illustrate  use  of  the  first  group.  Program  BGND3  (H20)  computes  vibrational 
populations  and  total  unit  volune  band  radiances  at  2  km  intervals  between  0  and 
150  km,  and  at  10  km  intervals  between  150  and  700  km.  Program  BGND9  integrates 
radiance  values,  mean  temperature,  and  total  number  of  molecules  per  square 
centimeter  in  each  vibrational  level  along  the  line-of-sight  for  two  endo- 
atmospheric  viewing  levels:  vertical  (0°)  and  horizontal  (90°)  and  for 
exo-atinospheric  limb  viewing.  The  computations  include  modifications  necessary 
tor  the  effects  of  optical  thickness  along  the  line-of-sight.  Program  SPCTRA 
computes  spectral  radiances,  corresponding  to  the  total  band  radiances  outputted 
oy  t*JND3  (H20)  and  BGND9,  between  250  and  6000  wavenumbers.  The  last  program 
i'KYPhN  plots  the  output  of  SPCTRA  for  a  maximum  of  five  altitudes,  one  at  a 
tune. 

Tne  following  is  a  user  oriented  discussion  and  description  concerning  the 
mechanism  and  input  data  required  to  make  a  successful  program  run.  A  detailed 
tlow  chart  of  each  program  is  also  included.  For  calculating  vibrational  pop¬ 
ulations  tne  flow  chart  for  the  H20  molecule  program  (BQ^D3)  is  used.  Only 
those  variables  v*iich  pertain  to  input  data  are  listed  and  derived.  A  complete 

list  of  the  variables  used  in  the  four  programs  and  their  related  subroutines  is 

r  3 1 

given  by  Degges  and  Smith1  . 

Subroutines  BNRAD2  (carbon  dioxide  program  BGND2)  and  BNRAD6  (ozone 
program  BGND6)  call  a  linear  equation  solver,  LEQTIF .  This  is  a  subroutine 
from  the  IMSL  library. 

This  library  is  licensed  for  use  at  AFGL  only,  and  we  cannot  supply  it 

with  the  proqrams.  The  IMSL  handbook  references  the  similar  subroutines 

(24  ] 

0EC0MP,  SOLVE,  and  SING  listed  and  described  by  Forsythe  and  Moler 
These  or  any  other  linear  equation  solvers  may  be  used  with  these  programs. 
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CALCULATE  VIBRATIONAL  POPULATIONS 


PROGRAM  BGND3 
(H20) 

The  main  program  BGND3  reads  and  checks  the  input  data  and  transfers 
control  to  subroutine  BNRAD3  which  computes  vibrational  population,  local 
optical  thickness,  and  volume  radiances  for  all  Included  bands  of  a  single 
molecular  species.  The  program  requires  as  input  the  temperature,  molecular 
concentration  data,  and  the  concentration  and  temperature  descriptors. 
PROGRAM 


Carbon 

Dioxide 

BGND2 

Water  Vapor 

BGND3 

Nitric 

Oxide 

BNGD4 

Ozone 

BNGD6 

Carbon 

Monoxide 

BNGD7 

See  Appendix  A  for  the  list  of  the  programs  and  the  data  on  the 
distribution  tape. 

PROGRAM  SPECIFI CATIONS 

There  is  a  different  program  to  be  used  for  each  molecular  species. 

PROGRAM  INPUT 

A.  Input  data  variables  or  constants  to  direct  program  flow  to  appropriate 
calculations  . 

Z FAC  Secant  of  solar  zenith  angle.  This  is  required  to  properly  treat 
absorption  of  incident  solar  radiative  fluxes . 

ICHEK  Used  to  print  diagnostic  inf orraat '  .  Normally  zero,  must  be  set 

to  l  (integer  one)  for  listing  of  optical  thicknesses  and  results 
of  Intermediate  radiative  transfer  calculations. 

ICOOL  Normally  zero,  set  to  integer  one  to  compute  and  print  atmospheric 
cooling  (or  heat ing)rates  for  each  infrared  band  included  in  the 
program . 

B.  Molecular  species  data:  Temperature,  Concentrations,  Band,  and  Level 
data.  See  Appendix  A.  Listed  below  are  the  variables  and  constants  in 
the  order  in  which  they  are  read  in. 

Integer  variable  used  to  indicate  conditions  of  computations  . 


NLTE 


“  1  -  Local  thermodynamic  equilibrium 

=  2  -  Collisional  excitation  only 

=  3  -  Night  conditions 

=*  4  -  Noon  conditions 

=  5  -  To  include  high  vibrational  temperatures  for  nitrogen 
T  Temperature  profile  (°K) 

Molecular  concentrations  (molecules/cm^ ) 

C0NN2  -  Nitrogen 
C0N02  -  Oxygen 
CONO  -  Atomic  Oxygen 
C0NCH4  -  Methane 
C0NC02  -  Carbon  Dioxide 
C0NH20  -  Water  Vapor 
CONNO  -  Nitric  Oxide 
C0NN20  -  Nitrous  Oxide 
CONO 3  -  Ozone 
CON CO  -  Carbon  Monoxide 

IGAS  Number  of  molecular  species  for  which  data 

is  to  be  read  in. 

XID  Array  of  molecular  species  names 

ICWC  Molecule  identifier  code 

1  -  Methane 

2  -  Carbon  Dioxide 

3  -  Water  Vapor 

4  -  Nitric  Oxide 

5  -  Nitrous  Oxide 

6  -  Ozone 

7  -  Carbon  Monoxide 
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ALTM 


TYPE 

NLEV 

NBAND 

AMAS 

DEXT 

DEXB 

LSC 

LCC 

RCM 

BV 

JGAS 
J  LEV  EL 


Maximum  altitude  of  the  molecular  species 
for  which  total  band  radiances  are  calculated  . 

Also  used  for  calculating  the  number  of 
altitude  levels.  The  values  In  (km)  are: 

120  -  Methane 

700  -  Carbon  Dioxide 

700  -  Water  Vapor 
700  -  Nitric  Oxide 
150  -  Nitrous  Oxide 
700  -  Ozone 
700  -  Carbon  Monoxide 

Temporary  location  for  species  names  (See  XID) 

Number  of  vibrational  levels  for  this  species 
Number  of  bands  for  this  species 
Molecular  mass  of  the  molecule  (g/mole) 

Rate  coefficient  at  300°K  for  vibrational  exchange 
of  Nitric  Oxide  with  Nitrogen  or  Oxygen  for  the  most 
significant  vibrational  level  (cm^/sec) 

Rate  coefficient  at  300°K  for  vibrational  de-excitation 
of  lowest  Nitric  Oxide  vibrational  level  (cm3/sec) 

Integer  array  containing  level  type  information 
Integer  array  containing  vibrational  level  coupling 
to  Nitrogen  or  Oxygen 

Array  containing  energy  level  of  each  band  (cm“^) 
Rotational  constant  (cm-^) 

See  ICHK 
Level  number 
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LBC  Integer  array  containing  information  on  which 

levels  produce  a  band 

LBU  Integer  array  containing  band  type  information 

SIR  Array  of  band  strength  for  each  band 

( cm” 1 ( atm-cm)” * ) 

SFLUX  Array  of  solar  flux  values  for  each  band 

TFLUX  Array  containing  fluxes  from  lower  altitudes 

for  each  band 

PROGRAM  OUTPUT 

Write  level  and  band  data  to  File  A 


SUBROUTINE  BNRAD3 


PROGRAM  SPECIFICATIONS 

Called  by  BGND3 


INPUT  DATA 


None  (read  in) 


OUTPUT  DATA 

Write  vibrational  populations,  unit 
optical  thickness,  and  volume  radiances 
to  file  3 


DECK  SET-UP 

1.  Attach  program  BGND3  . 

2.  Attach  required  subroutines  -  PR0GL1B. 

3  .  Attach  data  file  . 

4.  Catalog  output,  TAPE3,  TAPE4  . 
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PROGRAM  BGND3 


3JPART 


Program  //I 


|  PROGRAM  BGND9 [ 


PROGRAM  BCUD3 


I  Set  constants  for  Input  I 
I  and  Output  Files  I 

I  NR,  NP,  NW  | 


1  Read 

I  NLTE:  Computational  control  parameters.  | 
I  Temperature  and  Molecular  I 

I  concentration  data  I 


SUBROUTINE  BNRAD3 


DO  130  NZ  =  1,  NALT 
Initialize  Population  array  -  VIBPOP 
assuming  equilibrium 


DO  110  NL  =»  2,  NLEVEL 
Calc 

Vibrational  Population  of  each  level 
relative  to  the  ground  . 

Store  the  sum  of  Relative  Populations 
in  the  variable  -  SUM 


Calc 

Ground  State  Population 


DO  120  NL  !  2,  NLEVEL 
Calc 

Vibrational  Population  of  each  level 
by  multiplying  Total  Concentration 
by  Relative  Concentration  divided  by 
SUM 


DO  135  NL  *  2 ,  NLEVEL 
Initialize  Vibrational  Transfer  Variables 


DO  160  NZ  =  1,  NALT 
Calc 

Rate  at  which  a  single  molecule  in  a 
given  level  is  excited  or  de-excited 
by  collisions  at  each  altitude 


FIGURE  2A.  Detailed  Flow  Chart  for  Program  BGND3 
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FIGURE  2A.  Detailed  FLow  Chart  for  Program  BGND3  (cont'd) 
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PROGRAM  BGND9 


Program  BGND9  computes  the  effective  temperature,  and  the  radiance 

values  for  two  endo-atmospheric  viewing  levels  (horizontal,  vertical) 
and  one  exo-atmospheric  viewing  leveL  (limb).  In  the  current  version 

radiance  values  are  computed  at  2  km  intervals  of  tangent  height  up  to 
150  km  and  at  10  km  intervals  above  150  km.  The  integration  of  radiance 

values  for  each  viewing  level  begins  at  the  viewing  level  .  At  each  point 
on  the  line-of-sight ,  the  contribution  of  the  local  volume  emission  rate  to 

that  received  at  the  viewing  point  is  computed  . 

PROGRAM 

See  Appendix  A  for  program  requirements  . 

PROGRAM  MODIFICATIONS 

A.  If  the  viewing  angles  other  than  horizontal,  vertical, 
and  limb  are  desired,  they  must  be  defined  within  the 
program.  See  the  variables  ANGLE(K),  ZANGL,  and  DZANG 
in  the  program. 

B.  DO  830  Loop  controls  the  tangent  height  range,  currently 
set  for  60  to  250  km.  See  the  equation  below  the  program 
statement  DO  830  K=*3l,  86  to  determine  how  to  calculate 
the  indices  to  obtain  the  desired  tangent  height  range  . 

PROGRAM  INPUT 

Output  files  of  program  BGND3  .  No  other  input  is 
required  . 

PROGRAM  OUTPUT 

Write  integrated  radiances,  effective  temperatures  and 
total  number  of  molecules  to  file  3. 
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PROGRAM  BGND9 


FIGURE  3.  Detailed  Flow  Chart  for  Program  BGND9 
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PROGRAM  SPCTRA 


Program  SPCTRA  computes  spectral  radiances,  corresponding  to  total  band 
radiances  calculated  in  programs  BGND3  and  BGND9 ,  between  250  and  6000 
wavenumbers  .  The  program  in  its  present  form  can  calculate  spectral  radiance 

for  only  one  viewing  angle  at  a  time.  If  required  an  optional  filter  sub¬ 
routine  may  be  attached  to  the  program. 

PROGRAM 

See  Appendix  A  for  the  list  of  the  programs  on  the 
distribution  tape  . 

PROGRAM  MODIFICATIONS 
None 

PROGRAM  INPUT 

A.  Output  files  from  programs  BGND3  and  BGND9  . 

HTS  Altitude  (km) 

RAD  Total  band  radiances  . 

RAD( 1 , I)  exo-atmospheric  limb  viewing  level 

RAD(2,I)  endo-atmospheric  vertical  viewing  level 

RAD(3,I)  endo-atmospheric  horizontal  viewing  level 

THT  Dummy  altitude  variable 

AVTCMP  Array  in  which  effective  temperatures 

are  stored  (°K) 

AVTEMP(l.I)  Limb  viewing 
AVTEMP(2,I)  Vertical  viewing 
AVETMP(3,I)  Horizontal  viewing 


OOLCNT 


Array  In  which  column  counts  In  each 
vibrational  level  are  stored  . 


COLON T(  1,1)  Limb  viewing 
C0LCNK2,I)  Vertical  viewing 
00LCNT(3,I)  Horizontal  viewing 

B.  Input  data  variables  or  constants  to  direct  program  flow  to  appropriate 
calculations  . 


NLTE 

NALT 

JMAX 

KALI 

KMAX 

J  TMX 
K.LMX 
JLST 

JLDL 


See  PROGRAM  INPUT  list  of  variables  for 
program  BGND3 

Number  of  altitudes  for  which  results  are 
to  be  returned 

Number  of  angles  included  for  each  altitude 
Lowest  altitude  used  -  first  altitude  read  in 
Number  of  spectra  to  be  calculated  in  each 
trip  through  DO  Loop  500.  Five  is  maximum. 

Limit  of  DO  Loop  500  . 

Value  of  KMAX  on  the  last  trip  through  Loop  500  . 
Index  of  the  lowest  altitude  to  have  spectra 
computed  . 

Increment  of  altitude  index  . 
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KA,KB,KC, 
KD ,  KE 


Indices  for  selecting  the  desired 


viewing  level  . 

=  1  Limb 

-  2  Vertical  (0°) 

=  3  Horizontal  (90°) 

JLO  Dummy  variable.  Program  sets  JLST  to  JLO . 

JHI  Dummy  variable  .  Program  sets  JH1  to  JLO  +  1 

WTLO  Weighting  constant  =  1  .0 

WTHI  Weighting  constant  =  0  .0 

C.  SPECTRAL  LINE  DATA 

Spectral  line  data  is  read  in  subroutine  ROTATE. 
The  input  data  consists  of  line  positions  and 
strengths  for  vibrational-rotational  bands  of 
water,  carbon  dioxide,  ozone,  nitric  oxide,  carbon 
monoxide,  and  for  the  pure  rotational  spectrum 
of  water.  See  Appendix  B. 

PROGRAM  OUTPUT 

Write  spectral  radiances  for  each  altitude  to 
File  A  . 


DECK  SET-UP 


1  .  Attach  program  SPCTRA. 

2.  Attach  required  subroutines  -  PROGLIB. 

3.  Attach  data  file,  output  from  BGND3  and  BGND9  . 

4.  Attach  line  file. 

3.  Attach  input  data  file. 

6.  Catalog  output,  TAPE7  . 


PROGRAM  SPCTRA 


FIGURE  4 


.  Detailed  Flow  Oiart  for  Program  SPCTRA  (cont'd) 
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PROGRAM  TRY PEN 


The  plotting  program,  TRYPEN,  was  developed  under  a  previous  contract 
and  consists  of  a  main  program  and  six  subroutines,  LINSCA,  LOGSCA,  LINAXS, 
LOG  AXS ,  GRID,  and  NULINE . 

PROGRAM 

See  Appendix  A  for  program  requirements  . 

PROGRAM  MODIFICATIONS 
None 

PROGRAM  INPUT 

A.  Output  from  program  SPCTRA. 

YT  Spectral  radiances 

B.  Input  data  variables  or  constants  to  direct  program  flow  to  the 
appropriate  calculations. 

JMAX  Number  of  plots  per  run  . 

Eleven  is  maximum . 

IDENT  Job  name  or  job  identification. 

SX  Length  of  X-axis,  inches. 

SY  Length  of  Y-axis,  inches 

KALT  Lowest  altitude  plotted  . 

KDEL  Altitude  increment  between  plotted  curves. 

LL=0  Logarithmic  ordinate.  LCYCLE  will  be  the  maximum 

number  of  cycles  on  the  ordinate  scale  . 

LL=1  Linear  ordinate  scale.  LCYCLE=1,  forces  scaling 

of  all  sets  of  data  at  once.  LCYCLE=0,  scale  each 


set  of  data  separately.  Abscissa  scale  is  linear 


LCYCLE  Number  of  cycles  Co  be  plotted  on  ordinate. 

>0  No  more  than  LCYCLE  log  cycles  will  be  produced. 

=0  Input  value  of  MIN  will  be  used,  unless  it  is 

greater  than  or  equal  to  computed  value  of  MAX. 

<0  Exactly  IABS  (LCYCLE)  log  cycles  will  be  used. 

LW=0  Abscissa  -  Wavenumber. 

LW=1  Abscissa  -  Wavelength  (micrometers) 

LPHO"D=0  Ordinate  -  Power  (watts) 

LPHOT^l  Ordinate  -  Quanta  (photons)  per  second 

WMIN  Lower  wavenumber  limit  to  plot . 

Minimum  is  250  .0  . 

WMAX  Upper  wavenumber  limit  to  plot  . 

Maximum  is  6000  .0  . 

DELM  Lower  extension  to  plot  range  . 

HELP  Upper  extension  to  plot  range  . 

ISMP  Greater  than  zero  plot  JMAX  curves  on  the  same  plot  area 

PFAC  Scaling  factor  for  all  plots  . 

NXL  Number  of  characters  in  XBCD 

XBCD  Abscissa  label. 

NYI.  Number  of  characters  in  YBCD  . 

Y BCD  Ordinate  label. 

C.  Attach  plot  routines. 

PROGRAM  OUTPUT 

Program  will  generate  JMAX  number  of  plots, 
plotting  one  curve  (altitude)  per  graph. 

After  JMAX  plots  have  been  plotted  and 
labeled,  the  program  will  then  generate 
a  final  plot  with  JMAX  number  of  curves 
on  one  graph. 
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DECK.  SET-UP 


I  .  Attach  program  TRYPEN  . 

2.  Attach  required  subroutines  -  PROGLIB. 

3 .  Attach  plotting  routines  . 

4.  Attach  data  file,  output  from  SPCTRA. 
3.  Attach  input  data  file. 


i 


PROGRAM  TRY  PEN 


Read 

Plotting  parameters  and 
program  control  constants 


Calc 

XLOW 

-  Lower  extent  of  plot  range  . 

XHIGH 

-  Upper  extent  of  plot  range  . 

INCR 

-  Array  index  of  lowest  frequency 
to  be  plotted 

NX 

-  Number  of  spectral  radiance 
data  points  to  be  plotted  per 
plot  (altitude)  . 

NY 

-  Total  number  of  spectral 
radiance  data  points  to  be 
plotted  for  JMAX  number  of 
plots  (altitudes)  . 

DO  110  1=1,  2876 


Calc 

Frequency  array  -  XT 
XT(1)  -  250  .0  CM-1 

X T( 2876)  =  6000  .0  CM-1 


DO  114  I  =  1,  NX 

Store  frequencies  to  be  plotted  -  X 
X(l)  -  Lowest  frequency  to  be  plotted 
X(NX)  -  Highest  frequency  to  be  plotted 


FIGURE  5.  Detailed  Flow  Chart  for  Program  TRYPEN 
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v. 


DO  120  J  =  1,  JMAX 
Calc 

Array  indices  to  store  JMAX 
number  of  spectral  radiance 
data  arrays  in  a  single 
subscripted  array  -  Y  . 


DO  118  I  -  1,  NX 
Store  spectral  radiance  data 
in  Y  arra 


LPHOT 


Convert  spectral  radiance  -  P  (Watts) 
to  Quanta  (Photons)  per  seconds _ 


If 

I  =  NX 


Calc 

Number  of  elements  in  Y  array 
NXP  =  NX  +  2 
X(NX  +  1)  =  XLOW 
X(NX  +  2)  =  XHIGH 


FIGURE  5.  Detailed  Flow  Chart  for  Program  TRYPEN  (cont'd) 


DO  260  J  =  1.  JMAX 


Scale  each  set 
of  data 


Scale  all  sets 
of  data  at  once  | 


CALL  LINAXS  -  Draw  abscissa  axis 
CALL  LINAXS  -  Draw  ordinate  axis 
CALL  NULINE  -  Plot  data 


Plot  JMAX  number  of  curves  (altitudes) 
on  one  graph 

CALL  LINAXS  -  Draw  abscissa  axis 
CALL  LINAXS  -  Draw  ordinate  axis 
CALL  NULINE  -  Plot  data 


Rewind  spectral  radiance 
data  file 


CALL  ENDPLT 


FIGURE  5.  Detailed  Flow  Chart  for  Program  TRYPEN  (cont'd) 
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4.  OZONE  RADIANCE  MODELING 

Previous  versions  of  the  High  Altitude  Infrared  Radince  Model  have  pro¬ 
vided  for  only  one  band  of  ozone,  the  fundamental  transition  centered  at 
9.60  micrometers.  This  band  was  assumed  to  be  the  sole  channel  for  chemi¬ 
luminescent  radiation,  and  the  model  calculations  agreed  with  rocket  measure¬ 
ments  of  total  radiation  in  the  9  to  11  micrometer  spectral  region  to  within  a 
factor  of  about  two.  It  has  been  known  for  some  time,  however,  that  the  ozone 
is  formed  in  excited  levels  of  the  and  other  modes  in  the  three  body  reac¬ 
tion  of  atomic  and  molecular  oxygen.  Radiation  from  these  higher  levels  falls 
at  longer  wavelengths  up  to  about  12  micrometers,  partially  filling  in  the 
"window"  region  between  the  strong  ozone  fundamental  radiance  at  9.6  micro¬ 
meters  and  the  CC^  bands  be ween  12  and  16  micrometers.  Provision  for  these 
"hot"  bands  has  been  added  to  the  ozone  module  and  the  method  of  calculation 
is  outlined  here.  While  some  aspects  are  given  in  geater  datail  than  has  been 
done  for  the  other  molecules,  the  sketch  below  is  incomplete  and  requires  re¬ 
ference  to  earlier  work,  particularly  Ref.  3. 

f  231 

We  have  relied  heavily  on  Green  et  al.'s1  analysis  of  rocket  probe 

investigations  of  atmospheric  infrared  emissions  in  developing  the  current 

ozone  radiance  model  in  program  BNGD6.  The  references  in  the  tables  to  PSI 

refer  to  this  report.  In  modeling  the  infrared  radiance  field,  we  have  used 

r  i9i 

data  from  McClatchey  et  all  and  subsequent  updates.  Line  strengths  and 
positions  are  not  available  for  mode  bands  4-3,  5-4,.... and  we  have  in¬ 
cluded  them  by  extrapolating  from  the  1-0,  2-1,  and  302  data,  assuming  con¬ 

stant  dipole  moment,  but  correcting  for  altered  wavelength.  Table  11  lists 
the  vibrational  levels  included,  and  Table  12  lists  bands  used  in  the  model. 
Note  that  while  all  are  used  in  computing  vibrational  populations,  not  all  are 
available  for  spectral  radiance  calculations.  The  computational  scheme  is 
outlined  in  equations  (1)  -  (17).  We  use  square  brackets  []  to  indicate 
concentrations,  molecules  era-'*.  Parentheses  are  used  to  enclose  vibrational 
state  assignments.  (Oil)  indicates  no  excitation  of  the  mode,  and  one 

vibrational  quantum  in  each  of  the  V2  and  modes. 
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03(k)lk))J 


{  [03(000)]  (TF000010  +  CR000010) 

[03 (001)  ]  CR001010  +  [O3U00)]  CR100010 

+  [O3 (020) ]  (A020010  +  CR020010^ 

+  [03 (011) j  A011010  +  [O3(110)]  A110010 

+  [03(111)]  A111010}/(A0I0000  +  CR0i0000 
+  PR  4*  PR  4*  PR  4-  TR 

010001  010100  010020  010020 

4>  TF  4-  TF  +  TP  \ 

010011  010110  010111' 


Here  Tlr’000010  is  tne  rate  at  which  a  molecule  in  the  ground  state  (000)  is 
excited  to  the  mode  (010)  by  absorption  of  radiation  from  the  earth's  sur¬ 
face,  from  radiation  emitted  by  other  ozone  molecules  in  the  010-000  transition, 
and  (daytime  only)  from  the  sun.  The  collisional  excitation  rates  CR^  and 
guenching  rates  CH^  are  given  in  Table  13. 


[03  00J)j 


{  [O3(000)j  (TF000001  +  CR000001  +  cx000001> 
+  [O3(010)l  CR010(J01  +  IO3(100)]  CR100001 

+  “VM2>1  <*002001  +  CR002001)  +  (O3(101)1 
+  [Oj (111)1  A1U003  +  [0^  (003)  1  A003001 

»  ‘(o!  +a001000 


XJXJ  X 

’’’  CK001000  +  CXi 
+  CR^i,««  +  TF, 


t 

ana 


The  new  symbol  CF  introduced  here  is  the  rate  of  formation  of  ozone  in 
001 

tne  00I  state,  asstmed  to  be 

2.32  X  10"35  exp(510/T)(  [N2]  +  [02]>[021[0] 

The  assuned  rates  for  chemical  formation  in  this  and  other  vibrational 
states  are  presented  in  Table  14. 
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[Oj (100)1  = 


[0^020)]  = 


[O31011)J  = 


lo3ni0)]  = 


[U^ (002) j  = 


[U3(101)J  = 


[O3(000)]  (TF000100  +  CR000100) 

[°3(010)]  CR010100  +  [O3(001)]  CR00n00 

[03 (101) J  (A101100  +  CR101100) 

[03(200)]  A  +  f03 (HD]  Anil00> 
(1.0  X  10  [0]  +  A100000  +  CR100000 

CR100010  +  CR100001  +  TF100101 
CR100101  +  TF100200  +  TF100111) 

[°3(0i0)]  (tf010020  ♦  cr010020) 

(yX10  [O]  +  A020010  +  CR020010) 


(3) 


(4) 


1  t°3<000)]  TF000011  +  t°3(010)]  (TF0i00ii  +  CR0100H) 
+  [O3(001)]  CR0010U  +  [03 (111)]  Ani011) 

/  (1.0  X  10  [0]  +  A011000  +  A011010  +  CR011010 

+  CR011111  +  TF011111) 


(5) 


F000110  +  [O3(010)1  (TF010H0  +  CR010110} 

:R100110  +  l°3 (IH)I  Amii0} 

•JT01  +  A  +A  +  PR 

L  J  110000  110010  110010 

TF  ) 

H0iir 


(6) 


1  CF002  +  [O3<000))  FF000002 
+  [O3(001)J  (TF001002  +  CR001002)  + 

+  [03 (003) ]  (A003002  +  CR003002^ } 

/  (A002000  +  A002001  +  CR002001  \?R002101 
+  TF002003  +  CR002003  +  1,0  X  10 


(O3(101) 1CR101002 


(7) 


+  (0. 


33(001)] (TF001101  +  CR001101) 

+  CR100101)  +  tO3(002)]  CR0021 01 

4*  A  +  PR 

L01000  101001  101001 

'  CR101002) 


(8) 
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10jU“»J  -  1  ‘O3l000>  TF000200  +  IO3(100)1(T!t00200  +CR100200» 
+  [O3 (101)  J  CRi0i200}  ^  ^*0  X  10  ^0]  +  ^00000 

+  A200100  +  CR200100  +  CR200101) 


[0^(111) j  =  i  [o3(000)  tf0001U+  Co3(010)]  tf0101U+  to3(00i)]  tf001U1 

+  [CV100)]  TF100111  +  tO3<011)]{TF0iim  +  CR0imi> 

+  [O3(110)]^(TF110ill  +  CR110111)  +  [O3(101)]  CR101111)> 

/  (1.0  X  10  [0]  +  Alil000  +  A111010  +  A111001  +  Ami00 

+  A111011  +  CR111011  +  A111110  +  CR111110  +  CR111101) 


[O3(003) J 


10^(004) J 


1  CF003  +  t°3 (000) ]  TF000003  +  [03(001)]  TF00i003 
+  [O3 (002) ]  (TF002003  +  CR002003^  +  [°3(004)]  (A004003+CR004003) } 

/  (1.0  x  10  [o]  +  a003000  +  a003001  +  a003002  +  cr003002 

+  CK003004)  ^11) 

1  CF004  +  [°3(003)]  (TF003004  +  CR003004} 

+  [G3(005)J  (A005004  +  OR005004) > 

/  (1.0  X  10-n[O]  +  A004003  +  TF004005 

+  CR004005  +  CR004005)  (12) 


[O3(005) j 


1  CF005  +  [°3  (004)  ]  (^004005  +  CR004005) 

+  [o3(006)]^(a006005  +  cr006005)  ) 

/  (1.0  x  10  [o]  +  a005004  +  tf005006 

+  CR  +  CR  \ 

005005  005006' 


[G3(006) J 


=  <  cf006  +  to3C00S)i  (tf005006  ♦  cr005006) 
+  [03 (007) J  (A007006  +  CR007006) } 

/  (1.0  x  10-  [0]  ♦  a006005  +  tf006007 
+  CK006005  +  CR006007) 


[O3(007)] 


1  Ct007  +  [°3 (006) ]  (TF006007  +  CR006007) 
+  [O3(008)]  (A00g007  +  CR008007)} 

/  (1.0  x  10  [o]  +  a007006  +  tf00?008 

+  CR  +  CR  \ 

007006  007008' 
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t03 (aas>  j 


=  1  CF008  +  t°3<007»  ‘TF007008  +  CR007008> 
+  [O3(009)J  (A009008  +  CR009008H 

/  (1.0  X  10"U[OJ  +  a008007  t  tf008009 

+  CR  4*  CR  \ 

008007  008009' 


[O3(B09)1 


1  CF009  +  tOo(008)]  (TF00800g  +  cr008009) 
/  (1.0  X  10  [0]  +  ^009003  +  CR009008 

+  CR008007  +  CR008009) 


(16) 


(17) 


The  excitation  and  quenching  cycles  are  incomplete  in  this  set  of  equations. 
Further,  the  temperature  dependence  of  the  collisional  processes  is  not  known. 

In  the  program  itself,  the  notation  is  different.  Instead  of  using  a  vib¬ 
rational  level  designation  such  as  011,  we  use  a  level  number  06.  Then  the 
subscripts  may  be  included  in  a  six  or  seven  character  name  for  each  rate  co¬ 


efficient.  Thus  the  vibrational  excitation  rate  CR 


002003 


is  given  the  FORTRAN 


name  CK0812  and  the  radiation  absorption  coefficeint  is  called  TF0812. 

002000 

hinstein  A  coefficients  are  carried  in  array  EA  and  are  numbered  in  the  band 


sequence  of  Table  2,  so  that  is  EA<24) • 
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TABLE  11.  VIBRATIONAL  LEVELS  INCLUDED  IN  OZONE  RADIANCE  MODEL 


LEVEL  NO.  STATE  DESIGNATION  ENERGY,  CM 

1  OOO  0.000 

2  010  700.932 

3  001  1042.084 

4  100  1003.141 

5  020  1399.275 

6  OH  1726.528 

1  HO  1796.261 

8  002  2057.892 

9  101  2110.791 

10  200  2201.157 

n  HI  2785.241 

12  003  3041.200 

13  004  3988.0 

14  005  4910.0 

15  006  5803.0 

16  007  6665.0 

17  008  7497.0 

18  009  8299.0 
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TABLE  12.  INFRARED  VIBRATION-ROTATION  BANDS  INCLUDED  IN  OZONE  RADIANCE  MODEL 


BAM) 

TRANSITION 

ABSORPTION 

EINSTEIN  A 

BAND  CENTER 

NO. 

DESIGNATED 

STRENGTH* 

COEFFICIENT 

WAVELENGTH, 

1 

010-000 

6 .283( -19 ) 

0.224 

14.27 

2 

001-000 

1 . 394( -17 ) 

11.00 

9.60 

3 

100-000 

6.71K-19) 

0.594 

9.07 

4 

020-010 

4 . 164{ -20) 

0.446 

14.32 

5 

011-000 

5 . 373{ -20) 

0.116 

5.79 

6 

011-010 

4  - 503( -19) 

10.39 

9.75 

7 

110-000 

2 . 266 ( -20 ) 

0.053 

5.57 

8 

110-010 

1 .104 (-20) 

0.291 

9.13 

9 

002-000 

1 . 110( -19) 

0.342 

4.86 

10 

002-001 

1 . 743( -19) 

20.72 

9.84 

11 

101-000 

1 . 134( -18) 

3.673 

4.74 

12 

101-001 

3  -  222 ( -21 ) 

0.437 

9.36 

13 

101-100 

6.1 52 ( -20) 

9.681 

9.92 

14 

200-000 

3.000( -20) 

0.106 

4.54 

15 

200- 100 

6.296(-21) 

1.176 

9.11 

16 

111-000 

2  - 320( -20) 

0.131 

3.59 

17 

111-010 

3 . 850{ -20) 

3.669 

4.80 

18 

111-001 

1 . 520( -22 ) 

0.053 

5.74 

19 

111-100 

2.650(-22) 

0.116 

5.94 

20 

111-011 

1.650(-22) 

0.593 

9.45 

21 

111-110 

2.500(-21 ) 

11.01 

10.11 

22 

003-000 

1.1 10( -19) 

0.746 

3.29 

23 

003-001 

2. 230 (-21) 

1.026 

5.00 

24 

003-002 

1 . 960( -21 ) 

30.44 

10.17 

25 

004-003 

2 . 120( -23 ) 

36.34 

10.50 

26 

005-004 

2.658( -25) 

43.08 

10.85 

27 

006-005 

3 . 609( -27 ) 

48.49 

11.20 

28 

007-006 

5.485( -29) 

52.71 

11.60 

29 

008-007 

9 . 495 ( -31 ) 

56.12 

12.02 

30 

009-008 

1 .872( -32) 

58.67 

12.47 

*Units  are  CM  ^/(MOLECULE-CM  at  296K.  At  present,  spectra  are  not  available 
for  bands  12,  15,  17,  18,  19,  20,  21,  23  and  24. 
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TABLE 

13.  COLLISIONAL  OE 

n 

OZONE 

RADIANCE  MODEL. 

ft 

(1) 

03(010)  +  M 

> 

ft 

CR010000  = 

y* 

1 

(2) 

03(010)  +  M 

CR000010  = 

pi  ... 


PSI(14) 


=  3.0  X  10'14{[N2]  +  [023  +  [0]) 


CR000010  =  exp(-700.932  C2/T)  CR010000 

a  i  j  Aftinr  i j  /  /  ^  \  nil 


PS  I  ( 7 ) 


C2  =  1.438786  K/(CM  ) ,  SECOND  RADIATI00N  CONSTANT 


:3) 


:4) 


03(001)  +  M  +  03(000)  +  M 

CR001000  =  2-9  X  10'15  T‘/?  ([N2]  *  C02]1 

03(000)  +  M  +  03(001)  +  M 

CR000001  =  exp( -1042.084  C2/T)  CR001000 


PSI(24) 


PSI (6a) 


(5) 

(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


03(001)  +  M  -*  03(010)  +  M 

CRqoioio  =  5-0  X  io'14([n2]  +  [02]  +  [0]) 

03(010)  +  M  +  03(001)  +  M 

^R010001  =  exp( -341  - 152  C2/T)  CRqqioio 

03( 100)  +  M  ♦  03(000)  +  M 
CR100000  5  X  10  ([Ng]  +  [02]) 

03(000)  +  M  -  03(100)  +  M 

^R000100  =  exp(-1103.141  C2/T)  CR^qqqqq 

03( 100)  +  M  ♦  03( 010)  +  M 

CRlOooio  =  5'°  X  10"14^N2^  +  +  [0]) 

03(010)  +  M  ♦  03( 100)  +  M 

CR010100  =  exP(-402-209  C2/T)  CR100010 

03(100)  +  M  +  03(001)  +  M 

CR10oooi  =  1’°  X  1°’U^N2]  +  C02]  +  C0]) 


PSI(12) 

PSI (16) 

PSI(  -  ) 

PSI(  -  ) 

PSI(13) 

PSI(15) 

PSI (9) 
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TABLE 

OZONE 

13.  COLLISIONAL  DE-EXCITATION  AND  EXCITATION  MECHANISMS 

RADIANCE  MODEL  (Continued). 

INCLUDED  IN  THE 

(12) 

03(001)  +  M  -  03( 100)  +  M 

CR001100  =  exP(“61-057  c2/t)  criooooi 

PS  I ( 8 ) 

(13) 

03 ( 020 )  +  M  ♦  03(010)  +  M 

CR020010  =  2CR010000 

PSI(  -  ) 

(14) 

03(010)  +  M  ♦  03(020)  +  M 

CR010020  =  exP(-698*343  ^/T)  CR020010 

PSI(  -  ) 

(15) 

03(0U)  +  M  >  03(010)  +  M 

CR011010  =  CR001000 

PSI(  -  ) 

(16) 

03(010)  +  M  ♦  03(011)  +  M 

CR010011  =  exP ( -1025 .596  C2/T)  CRQll010 

PSI(  -  ) 

(17) 

03(011)  +  M  *  03(001)  +  M 

CR011001  =  CR010000 

PSI(  -  / 

(18) 

03(001)  +  M  +  03(011)  +  M 

CR001011  =  exP(-684-444  C27t^  CR01 1001 

PSI(  -  ) 

(19) 

03( 110)  +  M  *  03(010)  +  M 

CR011010  =  CR100000 

PSI(  -  ) 

(20) 

03(010)  +  M  +  03(110)  +  M 

CRoiono  =  exp(-1095*329  CRnooio 

PSI{  -  ) 

(21) 

03 ( 110)  +  M  -  03( 100)  +  M 

CR1 10100  =  CR010000 

PSI(  -  ) 

(22) 

03( 100)  +  M  +  03( 110)  +  M 

CR100110  =  exP(-693*120  C27T^  CR110100 

PSI(  -  ) 
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TABLE  13.  COLLISIONS  DE-EXCITATION  AND  EXCITATION  MECHANISMS  INCLUDED  IN  THE 
OZONE  RADIANCE  MODEL  (Continued). 


(23)  O3 ( 002 )  +  M  ♦  03 ( 001 )  +  M  PSI(23) 

CR002001  =  2  CR001000 

(24)  03(001)  +  M  ♦  O3 ( 002 )  +  M  P$I(  -  ) 

CR001002  =  exp(-1015.808  C2/T)  CRq02001 

(25)  03( 101 )  +  M  ♦  03 ( 001 )  +  M  PSI(  -  ) 

CR101001  =  CR 100000 

(26)  03(001)  +  M  -  03( 101 )  +  M  PSI(  -  ) 

CR001101  =  exp(-1068.707  C2/T)  CR1(moo 

(27)  03( 101 )  +  M  ♦  03 ( 1 00 )  +  M  PSI(  -  ) 

CR101 100  =  CR001000 

(28)  03C 100)  +  M  -  03( 101 )  +  M  PSI(  -  ) 

CR100101  =  exp( -1007.650  C2/T)  CR10110o 

(29)  03( 101)  +  M  ♦  03(002 )  +  M  PSI(ll) 

CR101002  =  1,0  X  10"11(1:N2^  +  [02]  +  [0]) 

(30)  03(002 )  +  M  *  03( 101)  +  M  PSI(IO) 

CR002 101  =  exP^"52-899  C2/J)  CR101002 

(31)  03(200)  +  M  -  03( 100)  +  M  PSI(  -  ) 

CR200100  =  2  CR 100000 

(32)  03(100)  +  M  ♦  03(200)  +  M  PSI(  -  ) 

CR100200  =  exp( -1098 .016  C2/T)  CR2001Q0 
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TABLE  13.  COLLISIONAL  DE-EXCITATION  AND  EXCITATION  MECHANISMS  INCLUDED  IN  THE 
OZONE  RADIANCE  MODEL  (Continued). 


(33)  03(200)  +  M  ♦  03( 101 )  +  M  PSI(  -  ) 

CR200101  =  CR 100001 

(34)  03( 101)  +  M  ♦  03(200)  +  M  PSI(  -  ) 

CR101200  =  exP(~l898,888  ^2/T)  £^200101 

(35)  03( 111 )  +  M  ♦  03(011)  +  M  PSI(  -  ) 

CR111011  =  CR100000 

(36)  03(011)  +  M  ♦  03 ( 1 1 1 )  +  M  PSI(  -  ) 

cRoi  1111  =  exP("^878*7^8  c2/t)  CR111011 

(37)  03( 111)  +  M  *  03(110)  +  M  PSI(  -  ) 

CR111110  =  CR001000 

(38)  O3CIIO)  +  M  *  03( 111 )  +  M  PSI(  -  ) 

CR1 1011 1  =  exP("998*980  ^/T)  CR111110 

(39)  03( 111 )  +  M  ♦  03( 101 )  +  M  PSI(  -  ) 

CR111101  =  CR010000 

(40)  03(  101 )  +  M  ♦  03( 111 )  +  M  PSK  -  ) 

CR101111  =  exP^647-450  C2yT)  CR111101 

(41)  03(003)  +  M  -  03(002)  +  M  PSI(22) 

CR003002  =  3  CR001000 

(42)  03(002)  +  M  *  03(003)  +  M  PSI(  -  ) 

CR002003  =  exp(-983.308  C2/T)  CR003002 

(43)  03(000)  +  0  ♦  03(001)  +  0  PS  I (6b) 

CXC00001  ■  ‘-0  x  10-11 
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TABLE  13.  COLL ISIONAL  DE-EXCITATION  AND  EXCITATION  MECHANISMS  INCLUDED  IN  THE 
OZONE  RADIANCE  MODEL  (Continued). 


(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 


03(001)  +  0  -  03(000)  +  0 

CX001000  =  exp( 1042.084  C2/T)  CXQ00001 

03(004)  +  M  -  03(003)  +  M 
CR004Q03  =  4  CR001000 

03(003)  +  M  *  03(004)  +  M 
CR003004  =  exp  *~986,8  C2/T^  CR004003 

03(005)  +  M  *  03(004)  +  M 
CR005004  =  5  CR001000 

03(004)  +  M  -  03(005 )  +  M 

CR003004  =  exp  f'922-0  C2/T)  CR005004 

03(006)  +  M  ♦  03( 005)  +  M 
CR006005  =  6  CR001000 

03 ( 005)  +  M  -*  03 ( 006 )  +  M 

CR005006  =  exp  (_893-°  C27T^  CR006005 

03(007)  +  M  *  03(006)  +  M 
CR007006  =  7  CR001000 

03(006)  +  M  -  03(007)  +  M 

CR006007  =  exp  ("862-°  C27T^  CR007006 

03 ( 008)  +  M  -  03(007)  +  M 
CR008007  =  8  CR001000 

03(007)  +  M  +  03(008)  +  M 
CR007008  =  exp  (_832-°  C2/T)  CR008007 


PSI(  -  ) 
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TABLE  13.  COLL IS ION AL  DE-EXCITATION  AND  EXCITATION  MECHANISMS  INCLUDED  IN  THE 
OZONE  RADIANCE  MODEL  (Continued). 

(55)  03(009)  +  M  ♦  03 ( 008 )  +  M 
CR009008  =  9  CR001000 

(56)  03(008)  +  M  -  03 ( 009 )  +  M 

CR008009  =  exp  {'802*0  C2/T)  CR009008 


TABLE  14.  CHEMICAL  FORMATION  RATES  FOR  OZONE  v3  EXCITED  LEVELS 


(1)  0  +  02  +  M  -  03(001)  +  M  PSI(lb) 

CF001  =  2.32  X  10‘35  exp(510/T)[0][02]([N2]  +  [Og]) 

(2)  0  +  02  +  M  -  03(002)  +  M  PSI(lc) 

CFoq2  =  1.85  X  10"35  exp(510/T)[0][02]([N2l  +  [02]) 

(3)  0  +  02  +  M  >  03(003)  +  M  PSI(ld) 

CF003  =  1,42  X  10"35  exp(510/T)C03[02KCN2]  +  [02]) 

(4)  0  +  02  +  M  -  03(004)  +  M  PSI(le) 

CF004  =  1.06  X  10~35  exp(510/T)[0][02]([N2]  +  [Og]) 

(5)  0  +  02  +  M  -  0 3 ( 00 5 )  +  M  PSI(lf) 

CF005  =  7.59  X  10'36  exp(510/T)[0][02]([N2]  +  [Og]) 

(6)  0  +  02  +  M  -  03(006)  +  M  PSI(lg) 

CF006  =  4.73  X  10"36  exp(510/T)[0][02]([N2]  +  [Og]) 

(7)  0  +  Og  +  M  +  0  3 ( 00  7 )  +  M  PSI(lh) 

CF00?  =  2.53  X  10'36  exp(510/T)[0][02]([N23  +  [Og]) 

(8)  0  +  02  +  M  -  03(008)  +  M  PSI(li) 

CFoo8  =  8,80  x  10~37  exP< sio/t) C03C023( [n23  +  [o2]) 

(9)  0  +  Og  +  M  *  03(009)  +  M  PSI(lj) 


CFq09  =  2.20  X  10'37  exp(510/T)[0][02]([N2]  +  [Og]) 
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[Xiring  development  of  the  new  ozone  bands,  we  made  a  change  in  the  method 
or  iteration  in  computing  vibrational  populations.  Instead  of  using  excited 
state  concentrations  from  previous  cycles  to  compute  new  vibrational  popula¬ 
tions,  the  calculation  was  recast  as  a  set  of  simultaneous  equations  for  vibra¬ 
tional  populations  at  a  given  altitude,  and  all  excited  state  concentrations 
at  a  single  altitude  are  computed  simultaneously.  This  avoids  the  problem 


noted  earlier  for  CC^,  that  the  closely  coupled  levels  may  hold  each  other 


close  to  the  initial  conditions  assumed  and  not  reach  convergence  in  a  reason- 
aole  number  of  iterations.  We  have  retained  the  earlier  method  as  comment 
stataments  m  the  program,  to  porvide  a  link  between  the  equations  above  and 
tne  present  coding  of  the  model. 


APPENDIX  A 


The  programs  in  a  FORTRAN  V  format  and  the  data  for  calculating 
vibrational  polulations  are  stored  on  one  magnetic  tape.  The  tape  s 
nine  track  1600  CPI,  the  information  is  stored  as  ASCII  Code  80  characters 
per  record.  There  are  60  files  on  tape.  The  following  is  the  order 
in  which  the  programs  and  data  are  stored. 


File 

Program/Data 

Destription 

1 

BGND2 

Carbon  Dioxide  Program 

2 

BGND3 

Water  Vapor  Program 

3 

BGND4 

Nitric  Oxide  Program 

4 

BGND6 

Ozone  Program 

5 

BGND7 

Carbon  Monoxide  Program 

6 

BGND9 

See  Page  56 

7 

SPCTRA 

See  Page  59 

8 

TRY  PEN 

See  Page  66 

9 

PROGLIB 

Program  Subroutines 

10 

BKNEW 

Molecular  species  data  common 
to  all  vibrational  population 
programs 

For  the  25  model  atmospheres  listed  below  there  are  night  and 
daytime  conditions  for  each  model  .  The  first  model  listed  is  the  1976 
standard  atmosphere  .  The  remaining  24  models  are  defined  by  three  exo- 
atmospheric  temperatures:  600°K,  1000°K,  and  1500°K  .  Thus  the  1976 
standard  model  requires  two  files  of  data  storage  while  the  other  24  models 


require  six 

files  of 

data  storage  per  model  . 

FILE 

MODEL 

EXO- ATMOSPHERIC 

TEMPERATURE 

11,12 

1976 

Standard 

1000 

Atmospheric 

13,14 

15° 

Annual 

600 

15,16 

15° 

Annual 

1000 

17,18 

15° 

Annual 

1500 

19,20 

30° 

Summer 

600 

21  ,22 

30° 

Summer 

1000 

23,24 

30° 

Summer 

1500 

25,26 

30° 

Winter 

600 

27,28 

30° 

Winter 

1000 

29,30 

30° 

Winter 

1500 

31,32 

45° 

Spring/Fall 

600 

33,34 

45° 

Spring/Fall 

1000 

35,36 

45° 

Spring/Fall 

1500 

37,38 

45° 

Summer 

600 

39,40 

45° 

Summer 

1000 

41  ,42 

45° 

Sumrae  r 

1500 

43,44 

45° 

Winter 

600 

45,46 

45° 

Winter 

1000 

«# 

4> 

00 

45° 

Winter 

1500 

49,50 

60° 

Summe  r 

600 

51,52 

60° 

Summer 

1000 

53,54 

60° 

Summe  r 

1500 

55,56 

60° 

Winter 

600 

57,58 

60° 

Winter 

1000 

59,60 

60° 

Winter 

1500 

The  file  positions  are  listed  11,12,  etc.  The  first  file  number 
is  for  the  night  conditions  and  the  second  file  number  is  for 
the  daytime  conditions 


K9 


NOTE  : 


APPENDIX  B 


The 

spectral 

line 

data  tape  "LINE 

TAPE"  contains 

five  files  of 

data  . 

FILE 

1 

Carbon  Dioxide 

2 

Water  Vapor 

3 

Nitric  Oxide 

4 

Ozone 

5 

Carbon  Monoxide 

In 

each  fiLe 

the 

data  is  ordered 

into  bands  with 

increasing  frequency 

within 

each  band 

• 

Carbon 

Dioxide 

Band  Nc 

i .  No  .  Lines 

Band  No  . 

No  .  Lines 

Band  No .  No 

.  Lines 

1 

L53 

11 

113 

21 

99 

2 

136 

12 

235 

22 

92 

3 

276 

13 

125 

23 

83 

4 

137 

14 

no 

24 

81 

5 

121 

13 

76 

25 

39 

6 

128 

16 

75 

26 

81 

7 

229 

17 

278 

27 

79 

8 

104 

18 

99 

28 

85 

9 

247 

19 

93 

29 

81 

10 

127 

20 

251 

90 


APPENDIX  B  (coat’d) 


Rotational  Band  -  L 58 1  Lines 


Band  No  . 

No  .  Lines 

Band  No  . 

No  .  Lines 

Band  No  . 

No  .  Lines 

L 

1773 

6 

34  5 

11 

* 

2 

1057 

7 

333 

12 

* 

3 

1315 

8 

662 

13 

989 

4 

1727 

9 

440 

14 

1233 

5 

559 

10 

737 

*  Bands 

11  and  12  are  not 

on  tape  . 

Nitric  Oxide 

Band  No 

.  No .  Lines 

Band  No  . 

No.  Lines 

Band  No . 

No.  Lines 

1 

234 

2 

193 

3 

160 

Ozone 

Band  No 

.  No  .  Lines 

Band  No  . 

No.  Lines 

Band  No  . 

No  .  Lines 

1 

6340 

11 

2165 

21 

* 

} 

5813 

12 

* 

22 

1575 

3 

5203 

13 

1185 

23 

* 

4 

4591 

14 

1530 

24 

3518 

5 

1709 

15 

■k 

25 

1534 

6 

1544 

16 

1031 

26 

1534 

7 

2137 

17 

* 

27 

1534 

8 

901 

18 

* 

28 

1534 

9 

2164 

19 

* 

29 

1534 

10 

1534 

20 

* 

30 

1534 

*  Ba  nd 

is  not  on  tape 

O.i  i  6on 

Monox 1  do 

Band  No .  No.  Lines 

Band  No 

•  No  .  Lines 

Band  No 

.  No.  Lines 

2  56  3  26 
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